Existence of Long-Range Order for Trapped Interacting Bosons 
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We derive an inequality governing "long range" order for a localized Bose-condensed state, relating 
the condensate fraction at a given temperature with effective curvature radius of the condensate and 
total particle number. For the specific example of a one-dimensional, harmonically trapped dilute 
Bose condensate, it is shown that the inequality gives an explicit upper bound for the Thomas-Fermi 
condensate size which may be tested in current experiments. 
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The classic Hohenberg theorem Q employs the fact 
that in the thermodynamic limit of an infinite system of 
interacting particles of bare mass m, the relation (Ti = 
fc B = l) 



mT no 



(1) 



for the occupation numbers Nk =< bjjik, > of the plane 
wave states enumerated by k 0) holds (the Bogoliubov 
1/fc 2 theorem 0). The angled brackets here and in what 
follows indicate a thermal ensemble (quasi-) average [Q; 
no and n are the condensate density and total density, 
respectively. The relation ([j]) entails that in one and two 
dimensions, a macroscopic occupation of a single state, 
the condensate, is impossible: The inequality leads to a 
contradiction in dimension D < 2 due to the (infrared) 
divergence of the wave vector space integral of (Q) , which 
determines the number density of non-condensate atoms. 
Physically, long range thermal fluctuations destroy the 
coherence expressed by the existence of the condensate. 

The question of the applicability of the Hohenberg the- 
orem to systems displaying long range order has renewed 
interest for trapped Bose-condensed vapors of reduced 
dimensionality, which can now be realized using various 
techniques J|, |J, |[ |[ 0. In these systems, it is possi- 
ble to investigate in a controlled manner the influence 
of finite extension in different spatial directions on the 
existence of a condensate, i.e. a macroscopically occu- 
pied single state created by b Q . For the non-interacting 
case, it is not difficult to show that there can exist a 
macroscopic occupation of one or several states in a Bose- 
condensed, trapped vapor of any dimension, as long as 
particle numbers remain finite, simply by evaluating the 
Bose-Einstein sums for the occupation numbers (see, e.g., 
H). The difficulty comes in when interaction is turned 
on. An interacting gas has a behavior increasingly dif- 
ferent from the ideal gas the lower the dimension of the 
system || [lCj. Furthermore, the classification of excited 
non-condensate states by plane waves is not the suitable 
one in a trapped gas: The condensation, in the limit of 
large total particle number N , takes place primarily into 
a single particle state in co-ordinate space jnj, and con- 
densate and total densities have (in principle arbitrary) 



spatial dependence, no — > no(r), n — » n(r). While a 
recent proof by Lieb and Seiringer shows that ground 
state condensation is, in the thermodynamic limit, 100% 
into the state that minimizes the Gross-Pitaevskii en- 
ergy functional in three as well as in two dimensions p2| , 
the question of the existence of a condensate for general 
energy functionals, at finite temperatures and, in partic- 
ular, in the generic one-dimensional (ID) case remains 
open. 

In what follows, we will take account of the phe- 
nomenon of spatially localized Bose-Einstein condensa- 
tion, by deriving an inequality analogous to the integral 
of (Q), with no explicit dependence on any Hamiltonian 
which has velocity independent interaction and trapping 
potentials. At a given finite temperature, the maximally 
allowed condensate fraction is related to an effective cur- 
vature radius of the condensate and the total particle 
number. The inequality thus allows for concrete state- 
ments on the limiting size of quantum coherent systems 
of reduced dimensionality, and their rcalizability for given 
condensate and total density distributions. As an appli- 
cation to a specific example, we consider a ID harmoni- 
cally trapped, dilute Bose-Einstein condensate, and it is 
shown that the relation leads to bounds on the parame- 
ters of a Thomas-Fermi cloud which may be verified in 
present experiments. 

Our analysis is based upon the following decomposition 
of field and density operators into condensate and non- 
condensate parts: 

$(r) = $n(r)b a + 6$(r) 
p(r) = \<S> (r)\ 2 bib + 5p(r) 
8p(r) = <P* (r)bl5$(r) + h.c. + 8&(r)S$(r), (2) 

where $o( r ) is the single particle condensate wave func- 
tion (normalized to unity) and iVo| ( I > o( T *)| 2 is the conden- 
sate density. The commutation relation 



6${r),6&(r') = 6{r - r 1 ) - $ (r)*5(r'), (3) 



is obtained from the canonical commutation relations for 
the quantum fields $(r), $t(r). As a consequence of this 



2 



commutator, the mixed response (the "anomalous" com- 
mutator) is taking the nonlocal form 



W $ (r) [|*S(r)* (r') - 8(r - r')} . (4) 



Here, 
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carrying out the commutator, < b >= 
> has been used, where this assignment is 



valid to 0(l/yQVo). The first term in the brackets on the 
RHS is due to the second term in the commutator (||) . It 
vanishes if one takes the Bogoliubov prescription that bo 
and 6q be replaced by c-numbers. It is a well-established 
fact that this (standard) Bogoliubov approach violates 
particle number conservation |l3| , and it becomes appar- 
ent below that neglecting the second term on the RHS of 
(|J) would lead to numerically strongly different predic- 
tions for the allowed condensate fraction. 

We now make use of this form of the mixed response 
in the Bogoliubov inequality O, which reads 
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(5) 



for any two operators A and C, where T is the tempera- 
ture. It is valid for any many-body quantum system, for 
which the indicated thermal (quasi-)averages are well- 
defined. We choose the operators in relation (||) to be 
the smeared excitation and total density operators 



A = / d D rf(r)5$(r), 
C = I d D rg(r)p(r), 



(6) 
(7) 



where f(r) and g(r) are complex regularization kernels. 

Next we derive a sum rule for the denominator 
on the RHS of (^), analogous to the /-sum rule 



J-^ dui u>S(k,uj) — Nk 2 /m for the dynamic structure 
factor 113], but in co-ordinate space. Using that 



p(r),H(r)\,p(r')\) = -i([V r • j(r),p(r')]) 
= ^- (S(r' - r) ($ t (r')A r $(r) + h.c.) 

-A r 5(r' - r) ($t(r')*(r) + h.c.) ) , (8) 



we obtain that the double commutator equals 



C,H 



d D r{-A r g(r))g*(r)n(r), 



(9) 

where n(r) =< p(r) >=< 3>'(r) <&(r) > . Here, we used 
the continuity equation for the current density operator 
j = [$tV$ - (V$ t )$]/2im, together with the Heisen- 
berg equation of motion for p: [p, H] — idtp — — iV r • j, 
this last relation being valid for a Hamiltonian with no 
explicit velocitydependence in the interaction and exter- 
nal potentials |H |. 

We normalize the kernel f(r) to unity, i.e., 
Jd D r\f{r)\ 2 = 1. The left hand side of the Bogoli- 
ubov inequality (JsJ) is then bounded from above, using 
the Cauchy-Schwarz inequality, as follows 



d D rld D 



7(r)/V) ({**(»•), ***(»•')}) 



= N-N + - 



(10) 



where the quantity a < 1 is given by a = 1 — 
Jd D r /d D r7(r)/*(r')$o(r)$5(r') , and where N - N = 
fd D r < Sp(r) >= Jd D r < 5& (r) <S$(r) > is the excited 
number of particles. 

Using the anomalous commutator (Q), the Bogoliubov 
inequality (||) thus may be written in the following form: 



N- N > -- 
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\Jd D rg(r)f(r)Mr) ~Jd D rJ d D r' g(r)f(r>)\<S> (r)\ 2 <S> (r>) \ 
fd D r(-A r g(r))g*(r)n(r) 

I 



(11) 



It is important to recognize that the above relation is 
explicitly independent of the form of the excitation spec- 
trum of the system, due to the relation (|J). In particu- 
lar, the strength of the interaction enters only implicitly 
in the form of the condensate wave function and total 
density distribution. This is in contrast to previous con- 
siderations on "long range" order in Bose-Einstein con- 
densates, employing correlation functions jHJ |l5|], where 



use was made of the excitation spectrum of the system, 
with Bogoliubov-type or WKB approximations. 

Now choose the kernels to have the particular form 



h(r) 



fi 1 ^ 2 exp[ifc • r] (re V ) 



(r i 2> ) 
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(12) 



where £Iq is the volume of the domain T>q in which $o 
has finite support. Hence \g\ 2 is also, automatically, nor- 
malized to unity, / d D r\g\ 2 = f T>g d D r\$ \ 2 = 1- The 



choice of g is motivated by requiring that it contains only 
information pertaining to (the spatial extension of) the 
condensate. The kernels are a set of functions in which 
we vary k such that the RHS of (jll]) is maximized and 
the inequality is assuming its strongest form. 



After introducing the kernels (12) into the Bogoliubov 
inequality (jll[), we have our primary result: 



1-7*7 



N - N > 



$ (fc) 



TN - *o(fc)/ Po d D r n(r)\Mr)\ 2 

2ih L d D r<f> (r) [-^A r %(r)] n{r) 



(13) 



The Fourier transforms of single particle conden- 
sate density and wave function are defined to be 
no(fc) = J d D r |<f>o(r)| 2 exp[ifc ■ r] and $o(fe) 
| d D r $o(f") exp[ifc ■ r]. 

We define the effective radius of curvature of the con- 
densate to be 



-1/2 



(14) 



According to the above formula, R c is obtained by weigh- 
ing a quantity proportional to the kinetic energy density 
of the condensate with n(r) / (N /SIq) , i.e., the local den- 
sity relative to the average density, and finally integrating 
over the domain T>q. Using this definition, relation Jl3| ) 
reads 



1 — F 

F 



>IH!^c(fe)- 5 l-(i-|$o(fe)| 2 M 



N 



x 2 



where the functional C(k) is given by 



C(k) = 



n (fe)-*o(fe) / d u rmr)\$ (r)\- 



(15) 



(16) 



and the condensate fraction F — Nq/N; the de Broglie 
thermal wavelength AdB = y/2w/mT. We stress that the 
value of C(k) is strongly reduced as a consequence of the 
commutation relatio n (|3j ), which causes the second term 
under the square in (|16|). 

The relation (JT^) implies that for temperatures close 
to zero, such that FR 2 /\\ B remains large (in which case 
the second term in (|TJ) is negligible), the approach of 
NC- 1 (1-F) =C- 1 (N-N ) ocT" to zero with a power 
law has to fulfill a < 1 for complete Bosc-Einstein con- 
densation into a localized single state &o(r) to be pos- 
sible. This statement holds for arbitrary strength and 
form of the interaction, in any spatial dimension. 

The above relations (|l3|) and (|l5|) are general. The 
strongest result, i.e., constraint on the system parame- 
ters we may expect for a ID system, in analogy to the 



original Hohcnberg theorem (the 1/fc infrared divergence 
of the integral of ([!]) in one dimension). To demonstrate 
the meaning of the relation (13) explicitly, we thus now 
proceed by considering the example of axially symmet- 
ric, harmonically trapped gases in one dimension, in the 
currently experimentally accessible Thomas-Fermi limit. 
Consider the Thomas-Fermi wave function 



7 2 
z/ Tp 



1/2 



(17) 



This mean-field form of <I>o is valid if the ID scatter- 
ing length fulfills the strong coupling condition n|aiD| 3> 
1, and the Thomas- Fermi parameters are Ztf = 
(3Ad4/| aiD |)i/3 ; n p = [(9/64)7V 2 |a 1D |/d4]i/3 g|. the 

quantities d±_ z = (mcjj_ z ) -1 / 2 are the harmonic oscil- 
lator lengths. In the quasi-3D scattering limit, which 
has transverse length scale d± ^S> a s , the ID scattering 
length is given by <2id = ~(d 2 L /2a s )[l — Ca s /d±} 1 with 
C = 1.4603 Jl6|. We neglect the difference between n(r) 



and N\<f> (r)\ , i.e., take the limit of both sides of ( pL 3|) to 
linear order in 1 — F. Evaluating the elementary integrals 
involved, we obtain that (x = kZ^-p) 



C(x) 



3 / SIM 



cos a; 



27tt 2 Ji(x) 



128 



(18) 



where J\ (x) is a Bessel function of the first kind. We see 
from Fig. [I] that the function C is strongly peaked at its 
global maximum k m ~ 3.7 Z^p (A m ~ 1.7 Ztf), where 
C(k m ) ~ 1.54 x 1(T 2 . Using that R c = 4Z T f/3\/2^, and 
neglecting the second term on the RHS of (|l5|) , we obtain 

cllj k kj-rn 



A. 



dB 



< 6.0 y/N - N . 



(19) 



We compare the above relation with the experiment on 
a ID 23 Na condensate in S, where the (relatively moder- 
ate) parameters were a s — 2.8 x 10~ 3 /^m, d z = 11.2 /im 
(lu z /2it = 3.5 Hz), d± = 1.15 /im, and N ~ 1.5 x 10 4 , 
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FIG. 1: The function C(k) on the RHS of 
dimensional Thomas-Fermi case, from Eq. (118 



in the one- 



which result in Ztf — 150 fim. Inequality ( |l9| ) becomes 
•^tf/Acib ^ 7.3 x 10 2 -\/l — F. Temperatures in this ex- 
periment have been of the order T ~ 0(100 nK) [fl7|| , 
which gives AdB — lpm and Ztp/X^b — 1-5 x 10 2 . The 
parameters of [|| are thus consistent with ([TJj]), provided 
F is not too close to unity jlqj . Note that jl9| ) would 
be inconsistent with the parameters in that experiment, 
were it not for the commutation relation ([}]), due to im- 
posing the canonical commutation relations for the total 
quantum field. Neglecting the second term on the RHS 
of ([}]) leads to a decrease of the RHS of (|l9|) by about one 
order of magnitude. In this sense, particle number con- 
servation, which is violated by the standard Bogoliubov 
prescription, is necessary for the condensate to exist |L9| . 

A further decrease of the aspect ratio ui z /u>± = d\jd\ 
down to values of order 10~ 3 and lower is achievable, 
e.g., in optical lattices @, so that the condition (|l9| ) 
on the trap parameters should be experimentally veri- 
fiable within present technology. When ( |l9| ) ceases to 
be fulfilled, the system has to enter into a new (non- 
Thomas-Fermi) state. For very low densities, a "Tonks 
gas" (ID gas of impenetrable particles) may be formed, 
which has no condensate [|o[ [n| ^FJ . Another possibility 
is a condensate consisting of (overlapping) phase coher- 
ent droplets M, resulting in a reduced value of R c . 

The primary result of the present investigation, in- 
equality (|l3|), holds under the generic conditions that the 
Bogoliubov inequality (|5|) is valid and that the potentials 
in the Hamiltonian are independent of particle velocities. 
It is, furthermore, to be emphasized that the application 
of ( |l3] ) is by no means limited to (real) ground state forms 
of <I>o. It is also possible to employ this relation to ex- 
amine existence conditions for excited state condensates 
with complex <E>o, like single vortices or vortex lattices. 
While ( |l3| ) certainly cannot guarantee the existence of a 
given {N , &o(r), n(r)} state, it can rule out models for 
trapped Bose-condensed gases in various spatial dimen- 
sions which are inconsistent with the Bogoliubov inequal- 
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